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Control design plays a major part in electrostatic linear comb actuator 
(ELCA) systems. Most existing control methods of ECLA mainly focus on 
traditional PID method. However, proportional—integral—derivative (PID) 
controller is only appropriate for the reference to be known as constant. In 
this article, the problem of establishing the dynamic model and designing the 


controller is investigated for ECLA working in the time-varying reference. 
The tracking effectiveness is constructed for ECLA systems based on the 
consideration of roots of the equivalent differential equation in tracking error 
model. Moreover, the convergence velocity is also tuned by changing the 
parameters in the presented control scheme. Simulation results guarantee the 
performance of the nonlinear control scheme as well as dynamic model to be 
established. 
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1. INTRODUCTION 

Micro electro-mechanical system (MEMS) is a remarkable device and popularly considered in several 
directions such as robots, vehicles, aerospace, industrial electronics, medicine, etc. Micro actuators can be 
known as one of the MEMS devices as well as describe an efficient task creating a force to move the micro 
elements and micro systems. Electrostatic linear comb actuator (ELCA) is a remarkable instrument structure 
using the working description via electrostatic force for establishing linear motion. The ELCA systems have 
many advantages including high effectiveness, low power requirement, and simple structure. They are able to 
explain the fact of easy movement in MEMS elements, such as resonator systems [1], gyroscopes [2] and micro 
systems of the conveyer or transportation systems [3], micro motors or gripper [4], [5]. Several work of ELCA 
have been recently developed with many appropriate results including new structures, working effectiveness, 
manufacturing improvement, smart materials. Particularly, the fabrication effectiveness and the ELCA’s 
movement stability guarantees the operating effectiveness of the elements. The characteristic of suspension 
beams with a crab-leg, a fold [6], [7] was developed to decrease the stiffness and change the deviation. A 
different direction on improving deviation was addressed [8], by decreasing the movement of out-of-plane 
from that generates the in-plane deviation. The work in [9] presented an actuator property after inserting 
more a module of a sensor as well as a sub-actuator for stability of control systems. The effect of side etching 
on deviation of ELCA with disability rent increase exemption DRIE method was pointed out in [10]. By 
computing the ELCA deviation in design objective [3], [4], [10]-[12], it is necessary to utilize the theoretical 
method to be developed by a static equilibrium equation and decreases a design time. 
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The fact is that these elements manipulate in statement of motion by using AC voltage such as 
traditional waves. Hence, under the effect of force and damping follows an efficient property in estimating 
the deviation of ELCA, which was considered in [13], [14] with the EAs to be developed for moving 
resonators, gyroscopes. Moreover, the motion dynamic model including an efficient mass and factor of 
conversion of air damping was not considered in these references. To achieve better behavior, the 
consideration of using control scheme is becoming a promising method [8]. However, the controller of 
ECLA have been considered with the traditional PID method. It should be noted that conventional PID 
controller does not satisfy the tracking effectiveness under the time-varying reference. Moreover, nonlinear 
controller and optimal control scheme have been usually proposed for different systems, such as robotic 
systems, power electronics [15]-[25]. Motivated by the above-mentioned discussions, this article studies the 
dynamic model of ECLA and the corresponding controller. To improve the tracking effectiveness under 
time-varying reference, the controller is carried out based on the consideration of roots of differential 
equations. 


2. PROBLEM STATEMENT 

Figure 1 describes the structure of ELCA including a moving element © and fixed element (@ and 
©). The moving element with moving comb fingers is pressured by four EB systems driving in y-direction 
based on the electrostatic force (EF), one terminal of each beam is kept on @. The configuration of movable and 
fixed comb fingers can be seen with the same shape. The work in [4] shows the ECLA geometrical parameters 
in Figure 1, where b is a beam width, h is the device layer thickness, and gais the gap. 
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Figure 1. Structure of ELCA 
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Under the influence of a voltage V between © and © of the ELCA, the total force in y-direction can 
be computed as (1). 


F= nhE£o y2 1 
pama (1) 
Where £o = 8.854 x 10-°pF/um and £= 1; V is the moving voltage. The stiffness in y-direction with four 
suspension beams system is described as (2) [4], [10]. 


_ 4Ehb? 
i 


Where E = 169GPa. 
Due to the application of a DC voltage for ELCA, the deviation yo of ELCA (i.e. displacement of the 
shuttle) is directly computed from the equation of force balance: 


(2) 


h 4Ehb? 
R — Fa = or: 2y = E-o 
It implies that: 
_ née (L E 2 
y= te (t)'v o 
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In (3) shows the quadratic relation between the deviation of the ELCA and the DC moving voltage. The 
objective of this article is to establish the dynamic model of ELCA for designing the controller with the 
reference of output y(t) depending on time. 


3. DYNAMIC MODEL OF ELCA AND CONTROL DESIGN 

Under the AC moving voltage, it implies the movement of ELCA’s shuttle to be similar with 
the movement to be described in Figure 2. Where K, C, M are the coefficients of conversion stiffness, 
conversion damping and conversion effective mass, respectively; F,(t) is an electrostatic moving force. 


Fe(t) 


Figure 2. 1-DOF dynamic model of the ELCA 


The dynamic model in y-direction (Figure 1) is represented as (4). 
My + Cy + Ky = R(t) (4) 


For the requirement of calculating exactly the deviation of the ELCA using dynamic model, the factors K, M, 
C are necessary to be accurately identified. Moreover, the factor K in y-direction is given in (2) and the 
appropriate mass of movable elements, including the elastic beams, is transformed into shuttle’s motion with 
the principle of the kinetic energy of mass M being equal to all movable parts’s kinetic energy. It follows 
that: 


. 2 
M =M, +% f" 2) dx (5) 


“L 40 OW 


moreover, the related equation between y(x) and y(L) is given through the relation of deviations: 


3Lx?-2x3 : 3Lx? -2x3 . 
y = Sy) > y) = FID) (6) 
according to (5) and (6), it leads: 
13M 52bL 
M=M,+4 se = phh +) (7) 


where p = 2330 kg/m? is the silicon mass density; A; is the total top side area of the shuttle. Based on [15], the 
coefficient C in (4) can be calculated: 


C=C, +C +C 


where 
C = u (+e) (8) 
C, = ps (9) 
Cs xut (10) 


combination of in (8), (9) and (10), we have: 


L (Abt 4 2Ase | AttAss | 16h 
at ea dor er ena (11) 


considering the driving voltage with the square wave (Figure 3). 
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Figure 3. Square wave of driving voltage 


The force function F. (t) is represented as: 


nhEE£o 


V? 0<t<T/2 
EE =} 6 Sry (12) 
0  T/2<t<T 


yı = e F(A, cos( wt) + A; sin(wt)) + - (13) 


Here: 6 = £, w = Jw? — B?; wn = 5 With the initial condition: y(0) = Oand y(0) = 0, it follows: 


BF 
Ay = -p 42 = Er 


E The general displacement formula of the shuttle in this case can be expressed: 


y= BC = (xan banca) i 


~~ 4Ego \b 


considering the driving voltage with the sine wave function (Figure 4). 


Figure 4. V = Vo(1 + sin(Qt)) 


NhEEQ 


By setting: Fy = 


Vē, the moving force is given: 
Jo 


Fas = Fo [= + 2 sin( 2t) — +cos(20¢)| (15) 


in this case, the deviation of dynamic model (4) is represented: 


3Fo 2Fo 


a — 2) si = 
oR (RUM [(K — MQ*) sin( 2t) — CAcos(Qt)] 
S n 2 ; 
NEOA [(K — 4M4 )cos(29Qt) + 2CQ sin( 22t)] 
+e FB, cos( wt) + Bz sin(wt)| (16) 


Vis = 


according to initial condition: y,,(0) = Oandy,,(0) = 0, we have: 


B, = 2c K-4M02 3 

17 "0 {cman + 2(K-4M22)2+8C202 =} 

B, = “of 2cn? _ _2A(K-M2?) }+ £ B, 
w ((K—-4MR2)2+4C?02? (K-MR2)2+C22 w 


for the mathematical model (4) established aboveMy + Cy + Ky = F,(t). The control objective is to achieve 
the input signalF,(t) such that y(t) tracks y”°f (t). Due to the difficulties in finding the precise parameters 
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M,C,K ; it is assumed that the parameters M, C, Kare unknown and bounded by known constants Mo, Co, Ko, 
respectively. Let’s define that e(t) = {y"°f (t) — y(t)}, according to (4): 


d2 d d2 ref d ref 
MES +C% + Ke = (M +C + ky — FO) 
that: 
2 
M3 +C&4+Ke = u(t) (17) 


2.,,7ref ref 
therefore, the proposed controller is chosen as F; (t) = (m 2 re +C ae + Ky" ) — u(t), where: 


de 
u(t) =D% (18) 
with D to be chosen as: 
2 
From (17) and (18), it implies that: M vo +(C-D) z +Ke=0 (19) 


hence, we only choose the parameter D such that the corresponding polynomial Mx? + (C — D)x + K = 0 
have all solutions being on the left side of complex plane. 


4. THE SIMULATION AND EVALUATION OF RESEARCH RESULTS 
The following simulations are performed to validate the tracking problem of output with the reference 
being the sinusoidal signal. The parameters of the proposed controller (18) are given as M = 1,C = 10, 


D=9,K =1. Two cases of sinusoidal reference signals are discussed as yı ref(t) = 2 sin (t+), 


Yoret(t) = 2 sin (t + z), The simulations will focus on the command tracking effect of the proposed controller 


(18). It is clearly seen that, although the references are depended on time but the tracking effectiveness is 
guaranteed (Figures 5 and 6). Moreover, the performance of tracking problem can be adjusted by choosing the 
parameters in (8) based on the roots of corresponding polynomial Mx? + (C — D)x +K = 0. 
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Figure 5. The result of the output with the reference being y,-¢(t) = 2 sin (t + z) 
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Response of Output 


Time (s) 


Figure 6. The result of the output with the reference being y,er(t) = 2 sin(t) 


5. CONCLUSION 

The article studies the dynamic model and control design of ELCA systems. Through the method of 
solving differential equation, the control scheme is implemented for ECLA in the presence of time varying 
reference. The tracking problem is analyzed through considering the solution of differential equation. By 
carrying out two simulation cases, it is concluded that the tracking problem is absolutely guaranteed with 
convergence rate to be depended on controller parameters. For future work, the problem of practical ECLA 
systems can be developed. 
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